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5y 1. Introduction 

Both the AGN and the nuclear starburst activities in galaxies which require inflow of material toward the centre either 
to fuel the central black hole or to cause rapid burst of nuclear star formation can be accomplished by tidal interactions 
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Abstract. We present Giant Meterwave Radio Telescope (GMRT) observations of the Hi 21 cm line and Very Large Array 
(VLA) observations of the OH 18 cm line from the Seyfert 2 galaxy Mrk 1. Hi emission is detected from both Mrk 1 and its 
companion NGC 45 1 . The H i emission morphology and the velocity field of Mrk I are disturbed. We speculate that the nuclear 
activities of Mrk 1 are triggered by tidal interactions. We estimate the Hi masses of Mrk 1 and NGC 451 to be 8.0(±0.6)xl0'* Mq 
and 1.3(±0.1) x 10' Mg respectively. We have also detected the Hi 21 cm line and the OH 18 cm line in absorption toward the 
nucleus of Mrk 1 at a blueshifted velocity with respect to its systemic velocity indicating an outflow of atomic and molecular 
gas. Two OH lines, at 1665 and 1667 MHz, are detected. Each of the profiles of the Hi and OH absorption consists of two 
components that are separated by ~ 125 km s"'. Gaussian fitting gave dispersions of ~ 44 km s"' for both the components of 
the H I absorption. The profile of the OH absorption is qualitatively similar to that of the H i absorption. Both components of the 
OH absorption are thermally excited. The peak optical depths of the two components of the H i absorption are (7.3 ± 0.4) x 10 " 
and (3.2 ± 0.4) x 10"^. The corresponding peak optical depths of the 1667 MHz OH absorption are (2.3 ± 0.3) x 10"^ and 
(1.1 ± 0.3) X 10"^. The higher velocity components of the Hi and OH (1667 MHz) absorption lines are blueshifted from 
the [O iii]/i5007, [Oi]/i6300, and the systemic velocity by ~ 100 km s"', but are consistent with the [Oii]/i3727 velocity. We 
explain these velocity discrepancies as due to shock ionization of a region which is pushed forward due to shocks in front of 
the radio nucleus thereby giving apparent blueshift to Hi, OH, and [On] velocities. The optical depth ratios thi/Tq^' of both 
the components of the Hi and OH absorption are ~ 3, indicating their origin in dense molecular clouds. Using OH/Ay values 
for the Galactic molecular clouds, we obtain 9 < Ay < 90 toward the line of sight of Mrk 1. 

Key words, galaxies: active - galaxies: interactions - galaxies: individual (Mrk 1, NGC 451) ~ galaxies: ISM 



( Hernquist & Mihos 1995 ). It is not clear, however, in the case of Seyfert galaxies whether nuclear activities in these low luminos- 



ity active galactic nuclei (AGN) are due to interactions as found in QSOs, radio galaxies, and BL Lacs (see De Robertis et al. 1998 
for a review on the subject). It is generally accepted that interactions leading to mergers (bound interactio ns) may play a sig 
nificant role in triggering nuclear activities than unbound or hyperbolic encounters ( fDe Robertis et al. 1998 ). Interactions can be 



effectively traced via Hi 21 cm line emission from galaxies as Hi disks often extend well beyond the optical radii of galaxies 
where the disks respond quickly to gravitational perturbations. H i emission studies may be particularly useful since most often 



the H I morphology provides evidences of interactions which are undetectable at optical wavelengths (e.g., Simkin et al. 1987). 

H I in absorption can trace kinematics and distribution of atomic gas near the centres of active galaxies on the size scales 
of their background radio sources. The advantage of absorption studies is that it can detect relatively small quantities of gas 



irrespective of the redshift of the object. Recently, Gallimore et al. (1999) have found H i rich absorbing disks on the scales of 
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Table 1. Global properties of Mrk 1 



Parameter 


Value 


Reference 


Right Ascension (J2000) 


oi'n6™OT:i 


1 


Declination (J2000) 


33°05'22" 


1 


Distance (Mpc) 


68 


2 


Hubble type 


SB 0/a 


1 


Seyfert type 


2 


1 


Inclination 


45° 


3 


Optical diameter (kpc) 


9.0x5.3 


1 


Connected blue magnitude BZ 

& T 


14.53 


4 


Total blue luminosity (Lq) 


1.1 X 10'" 




Total H I mass (Mq) 


8.0 X 10** 


5 


Hi mass to blue luminosity ratio (Mq/Lo) 


0.07 


5 


Total FIR luminosity (Lq) 


1.7 X 10'° 


6 


1.4 GHz radio luminosity (W Hz"') 


4.2 X 10-2 


5 


Spectral index (S oc y-") (2.7 GHz - 10.5 GHz) 


0.8 


7 


X-ray luminosity (erg s" ' ) 


< lo-" 


8 


Systemic velocity (km s"') 


4780±2 


9 


[Olll] A 5007 velocity (km s"') 


4822±25 


10 


[On] A 3727 velocity (km s"') 


4697±25 


10 


[Ol] A 6300 velocity (km s"') 


4817±25 


10 


Mean velocity of Hi emission (km s"') 


4780 


5 


H2O maser velocity (km s"') 


4868 ± 1 


3 


Mean velocity of CO emission (km s ') 


4850 


11 



Notes: Hq = 75 km s ' Mpc ' . The velocity definition is optical and Helio-centric. 
1: Markarian et al. (1989); 2: White et al. (1999); 3 Braatz et al. (1997); 4: NED (NASA Extra- 
galactic Database); 5:This paper; 6:IRAS faint source catalog, (1990); 7: Dickinson et al. (1976); 8: 
Fabbiano et al. (1992); 9: Keel (1996); 10: De Robertis & Shaw (1990) 11: Vila-Vilaro et al. (1998) 



a few hundred parsecs in several Seyfert galaxies. As a result of intense nuclear activities, gas in the central regions of active 
galaxies may be perturbed due to interactions of the radio p lasma with the surrounding ISM which may result in bulk outflows 
of material (e.g., Tadhunter et al. 2001 , Morganti et al. 1998 ). The molecular gas near the centres of active galaxies can be traced 
via 18 cm OH line in absorption. The 18 cm OH absorption line is sensitive to molecular gas in both the diffuse I SM and in the 
dark clouds with OH to H2 ratio being almost constant over a large variety of Galactic clouds ( Liszt & Lucas 1996 ). Studies have 
shown that chances of detecting OH absorption are higher in infrared luminous galaxies (Schmelz et al. 1986). 

In this paper, we present synthesis observations of the H 1 21 cm line obtained with the GMRT and the OH 18 cm line obtained 
with the VLA of the infrared luminous active galaxy Mrk 1 and its companion NGC 45 1 . The global properties of Mrk 1 are 
summarized in the next section. The details of observations and data analyses are given in Sect. 3. The results are presented in 
Sect. 4. Sect. 5 discusses the radio continuum properties, H i emission, and H i and OH absorption. The conclusions are in the 
last section. 



2. Global properties of Mrk 1 

The global properties of Mrk 1 are liste d in Table ^ Mrk 1 (NGC 449; B!^ = 14.53) is a member of a poor group (WBL 035) 
at a redshift of 0.017 ( |White et al. 1999| ). The other two members of this group viz., NGC 447 and NGC 451, are at projected 
separati ons from Mrk 1 of ~ 3 8 kpc and ~ 130 kpc respectively. Mrk 1 is classified as a SB O/a galaxy with a Seyfert type 2 
nucleus ( Markarian et al. 1989 ) with no signatures of interactions in the optical images. Mrk 1 is also a luminous Infrared galaxy 
(Lfir = 1.7 X IO^'Lq), indicating a high rate of star for mation. Mrk 1 is on e among 16 galaxies detected in the 22 GHz water 
megamaser em iss ion in a sample of 354 active galaxies ( Braatz et al. 1994 ). The nuclear optical spectrum of Mrk 1 studied by 
Koski (1978) and Weedman et al. (1968) shows strong emission lines typical of an active galaxy photo-ionized by ha rd contin- 
uum. The broa d lines indicative of a hidden Seyfert nucleus are not found either in the infrared ( [Veilleux et al. 1997) or in the 



polarized light (Kay et al. 1994). High dispersion spectroscopic observation of the [Oiii]/I5007 line by Bergeron & Durret (1987) 



shows a distinct blue asymmetry indicative of an outflow of gas. Keel (1996) suggested that the nuclear activities of Mrk 1 are 
due to an ongoing interaction with the nearby galaxy NGC 45 1 . 

The radio continuum emission from Mrk 1 is known to have a steep spectrum with a spectral index a (S oc v ") of 0.8 
(Dickinson et al. 1976). The 1.6 GHz EVN image (resolution ~ 30 pc) of Kukula et al. (1999) shows the nuclear emission to 
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Table 2. Observational parameters 



Parameter 


GMRT 


VLA 


Dates of Observations 


2000 Oct 25, 28 


2001 Mar 26, 27 


Pointing centre (RA J2000.0) 


01''16'"07:25 


0l''16"'07!25 


Pointing centre (Dec J2000.0) 


+33°05'22'.'2 


+33°05'22'.'2 


Observing duration (hrs) 


8 


4.5+ 


Range of baselines (km) 


0.1-25 


0.1-11 


Observing frequency (MHz) 


1395.00 


1640.22 


Bandwidth per IF (MHz) 


8.0 


6.25 


Number of spectral channels 


128 


128 


Polarizations 


2 


1 


Frequency resolution (kHz) 


62.5 


48.8 


Velocity resolution (km s"') 


13.7 


9.0 


Amplitude calibrator 


0137+331 


0137+331 


Phase calibrator 


0137+331 


0137+331 


Bandpass calibrator 


0137+331 


0137+331 



Usable time, see section 3.2 for details 



consists of an unresolved core surrounded by a weak diffuse emission with a total flux density of 34 mJy. The NVSS flux density 
at 1.4 GHz is 75.4 mJy 

The Arecibo observations by Hutchings (1989) detected Hi emission and blueshifted Hi absorption from Mrk 1. This single 
dish spectrum could not separate H i emission from Mrk 1 and NGC 45 1 . Observations with the Nobeyama Radio Telescope 
detected CO (J - 1-0) emission with a total flux integral of 11.5 + 1.6 K km s~' from the central 5 kpc region of Mrk 1 
( Vila-Vilaro et al. 1998 ). The search for the 1 8 cm OH absorption by Schmelz et al. (1986) with the Arecibo reflector resulted in 
a non-detection with an rms sensitivity to an optical depth of 0.02. 



3. Observations and data analyses 

3.1. The GMRT observations 

The GMRT observations of Mrk 1 were carried out in October, 2000. A summary of the main observational parameters are given 
in Table |[ At the time of the observations, the GMRT was not fully operational and hence not all 30 antennas were available for 
observations at any given time. Two runs of observations with 18-20 antennas, each with a field of view (FWHM) 24' centered 



on Mrk 1, were carried out on two different days. The GMRT has a mix of both short and long baselines (see [Swarup et al. 1991 
for more details), making it sensitive to diffuse emission of extent as much as 7' while having a maximum resolution of ~ 3" at 
1.4 GHz. The GMRT uses a 30-station FX correlator which produces complex visibilities over 128 spectral channels in each of 
the two polarizations. The bandwidth can be selected in multiples of 2 between 62.5 kHz and 16 MHz. These observations were 
carried out with a bandwidth of 8 MHz centered at 1395.0 MHz, which covered Hi velocities in the range 3730 - 5460 km s"^ 
with a velocity resolution of ~ 14 km s"'. 

The complex gains of the antennas were determined every 30 minute using observations of an unresolved nearby (~ 4.5°) 
source (3C 48) for 5 minutes. 3C 48 was also used for the flux and the bandpass calibrations. The data were reduced, following 
standard calibration and imaging methods, using the Astronomical Image Processing Software (AlPS) developed by the NRAO. 
The data were calibrated for the amplitude, phase, and frequency response for all antennas separately for each polarization. The 
flux density of 3C 48 was estimated to be 16.228 Jy at the observing frequency using the 1999.2 VLA flux densities of the 
standard VLA flux calibrators and the formula given in the AlPS task 'SETJY'. Due to the close proximity of 3C 48 to Mrk 1 
and based on some previous test experiments, we expect that the flux calibration is accurate to within 5%. 

A continuum data set was formed by averaging 80 line-free channels. The data were self calibrated in both phase and am- 
plitude. The resulting antenna gain corrections were applied to all channels. The continuum images were made using the self- 
calibrated averaged data from the line-free channels. The continuum flux density from each individual channel was subtracted 
in the (u,v) dataset by a Unear fit to the visibilities in the line-free channels. Since these observations were also sensitive to Hi 
emission, the data points were "natural-weighted" to enhance sensitivity to extended features. The resulting spectral cubes were 
CLEANed for signals greater than 4 times the rms noise in the channel images. The cube was blanked for emission below a level 
of 1 .5cr in the images after applying a Hanning smoothing of three velocity channels and Gaussian smoothing of five pixels (pixel 
size - 6") in the spatial co-ordinates. The zeroth and first order moment maps were generated from the blanked channel images 
containing H i emission and two additional channels on both the sides. 
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Fig. 1. Radio continuum images of the group WBL 035. The image in the left panel is at 1.4 GHz from the GMRT and that at the 
right panel is at 1.6 GHz from the VLA. The crosses mark the optical positions of the members of the group. The contours are 
drawn as -1, 1, 1.5, 2, 4, 8, 16, 32, 64, 128 in units of 1.1 mJy beam"' at 1.4 GHz and in units of 0.5 mJy beam"' at 1.64 GHz 
respectively. The synthesized beam, shown in the bottom left hand corner, is 30'.'75 x 25'.' 17, PA - +67.8°in the 1.4 GHz image; 
and 17'.'49 x 16'.'56, PA - -52.0°in the 1.6 GHz image. The images are corrected for respective primary beam attenuations. The 
rms noise is 0.35 mJy beam"' in the 1.4 GHz image and 0.19 mJy beam"' in the 1.6 GHz image. The peak flux densities in the 
images are 76.3 mJy beam"' at 1.4 GHz and 68.1 mJy beam"' at 1.6 GHz. 



3.2. The VLA observations 

The VLA 'B' configuration observations were carried out in March, 2001. The observational parameters are summarized in 
Table |[ The data were recorded in the lA correlator mode with a total bandwidth of 6.25 MHz divided into 128 channels. 
These observations covered a velocity range of 4400 - 5550 km s"' for the 1667 MHz line of OH with a velocity resolution of 
~ 9 km s"'. These observations also covered the 1665 MHz line of OH in the velocity range of 4050 - 5200 km s"'. A large 
fluctuation in the system temperature was noticed on the first day of observations which were carried out at a centre frequency 
of 1640.5 MHz. These fluctuations were later identified as due to strong signals from a external interfering source operating near 
the frequency of observations. The observations on the next day were carried out after reducing the front-end bandwidth from 
25 MHz to 12.5 MHz and shifting the centre of the band to 1640.22 MHz to avoid the external interference. The data from the 
first day of observations were discarded. 

The VLA observations and data analyses were carried out following the same strategy as adopted for the GMRT observations 
described in section 3.1. The flux density of 3C 48 was estimated to be 14.270 Jy at the observing frequency. The image cube 
was generated using "natural-weighted" continuum-free data to get maximum signal to noise ratio. The image cube was box-car 
smoothed along the frequency axis using a window of 3 channels and every second channel was discarded. The resulting image 
cube has a resolution of 5'.'43 x 4'.'93 x 27 km s"'. The continuum images were made using the self-calibrated averaged data 
from the line-free channels. 



4. Results 

4.1. Radio continuum 

The radio continuum images shown in Fig. [l| have a resolution of 30'.'75 x 25'.' 17 at 1.4 GHz, and 17'.'49 x 16'.'56 at 1.6 GHz. 
These images were made using only short (u,v) spacings to enhance the sensitivity to extended features. These images have 
an rms of 0.35 mJy beam"' at 1.4 GHz, and 0.19 mJy beam"' at 1.6 GHz. Continuum emission is detected from both Mrk 1 
and NGC 451. NGC 447 is marginaUy detected (~ 4cr) at both 1.4 and 1.6 GHz. The flux density of Mrk 1 is estimated to be 
76 + 4 mJy at 1.4 GHz, and 68 ± 3 mJy at 1.6 GHz. The spectral index between 1.4 and 1.6 GHz is 0.8. The flux density of 
NGC 451 is 3.3 + 0.5 mJy at 1.4 GHz and 1.8 + 0.4 mJy at 1.6 GHz. Mrk 1 remains unresolved down to a resolution of ~ 1 kpc. 



4.2. H I emission 

The H I cube was made with a resolution of 30'.'67 x 27'.'28 x 13.7 km s"'. The channel images have anrms of 0.92 mJy beam"'. 
The corresponding 3cr sensitivity in H i column density is 5.0 x 10'^ cm"^. The channel images of H i emission and absorption 
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Fig. 2. Channel images from GMRT showing H i column density contours in the velocity range 4807 km s ' - 501 1 km s ' . The 
crosses mark the optical positions of Mrk 1 and NGC 45 1 . Solid contours representing column density of H i emission are drawn 
at 3.6, 5.4, 7.2, 9.0, 10.8, and 12.7 xlO'^ cm"^. The negative contours (dashed curves) are drawn at 2, 3, 4, 5, 6 mJy beam"'. The 
HPBW of the synthesized beam (30'.'67 x 27 '.'28, PA = -80.4°) is indicated at the bottom left hand corner of the first channel 
image. The velocity resolution in the cube is ~ 13.7 km s"' . 



are shown in Figs. ^ and |] The results are summarized in Table || H i emission is detected in the velocity range of 4698 to 4984 
km s ' . Hi emission from Mrk 1 appears from 4698 to 4848 km s ' . The flux integral ( J 5 dV) of H i emission in this range is 
0.73 + 0.05 Jy km s"' which at the distance of Mrk 1 coiTesponds to a total Hi mass of 8.0(+0.6) x 10** Mq for Mrk 1. The 
estimated HI mass is a lower limit due to the effects of absorption. H i emission from NGC 45 1 is detected from 4807 to 4984 
km s"' with a total flux integral of 1.20 + 0.10 Jy km s"' corresponding to a Hi mass of 1.31(+0.11) x lO'^ M© assuming a 
distance to NGC 451 of 68 Mpc. The summed flux integral of Mrk 1 and NGC 451 is 1.93(+0.1 1) Jy km s ', which is consistent 
with the value obtained by the single dish observations of Hutchings (1989) See Fig. ^for a global H i profile. 
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Fig. 3. Channel images showing column density contours in the velocity range 4589 km s ' - 4793 km s The contour levels 
are the same as in Fig. |[ The H i absorption is seen toward Mrk 1 as dotted contours. 



The moment zero map shown in Fig. g indicates that the H i emission from Mrk 1 is distributed in three clumps with almost all 
H I seen outside the optical extent of Mrk 1 . The individual clumps having velocity dispersions of 30 to 60 km s ' are distributed 
over an extent of ~ 30 kpc. The velocity field of Mrk 1 shown in Fig. || indicates that there is a smooth rotation of H i from one 
end to the other. 

The Hi emission from NGC 451 shown in Fig. H looks like that of a disk galaxy with a total projected velocity width of 
170 km s The Hi diameter of NGC 451 is ~2Q kpc which is about twice that of the optical disk. The global parameters of 
NGC 451 given in Table|| were derived from a fit to the velocity field made using a higher resolution (17" x 14") Hi cube which 
is not shown here. 
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Fig. 4. Global Hi emission profile of Mrk 1 and NGC 451 from GMRT. The flux integral is 1.93 
consistent with the single dish observations of Hutchings (1989) . 
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Fig. 5. The column density contours of the total H i image from GMRT of Mrk 1 (top) and NGC 45 1 (bottom) overlaid upon 
the grey scale optical image from the DSS (blue). The contour levels are 0.3, 0.8, 1.3, 1.8, 2.3, 3, 4, 5, 6, 7, 8, and 9 in units 
of 10^" cm 2. The HPBW of the synthesized beam, shown in the bottom right hand corner, is 30'.'67 x 27'.'28, PA = -80.4°. 
Although the H i emission is surrounding Mrk 1, Hi absorption (marked as white circle) is detected toward Mrk 1 indicating the 
presence of cold H i gas in front of it. 



4.3. Hi absorption 

The channel images shown in Fig. |]also show H i absorption from Mrk 1 (dotted contours at the location of Mrk 1) in the velocity 
range 4589 km s"' to 4752 km s"' . The H i absorption spectrum shown in Fig. ^ is extracted at the radio position of Mrk 1 from 
a Hi cube made with a resolution of 6'.'23 x 4'.'39 x 13.7 km s"'. The rms in this cube was 0.6 mJy beam"'. The spectrum 
shows a broad multi-component absorption in between the velocities 4500 km s"' and 4800 km s"'. Two Gaussian components 
were fitted to the Hi absorption profile. The resulting parameters of the fit are given in Table Q The peak optical depths of the 
two components are 0.073 + 0.004 and 0.032 + 0.004 respectively and the velocity dispersions are ~ 44 km s ' for both the 
components. The column density of H i is estimated using the relation A^H I = 1-82 X 10'** x iTspi„/f) Jt dv cm ; where Tsnin 
is the spin temperature of H i in kelvin, / is the covering fraction of H i gas, J t dv is the velocity integrated optical depth in km 
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Fig. 6. The velocity fields of Mrk 1 and NGC 451 from GMRT are shown as contours and in grey scale. 



Table 3. Radio continuum and H i emission results 



Parameter 


Mrkl 


NGC 451 


srS.(mJy) 


76 ±4 


3.3 ±0.5 




68 ±3 


1.8 ±0.4 


^y(H I emission) (laji^-l) 


120 


170 


Systemic velocity (km s"') 


4780 ± 13"'" 


4897 ± 2* 


Maximum rot. velocity (km s"') 




140* 


Inclination (degree) 




31 ± 10* 


H I extent (kpc) 


~ 30 


~ 20 


/ 5 (Hi emission) dV(JykmS-') 


0.73 ± 0.05 


1.2 ±0.1 


Hi mass (lO** Mq) 


8.0 ±0.6 


13 ± 1.0 



t: Estimated from the mean velocity of Hi emission 

p. Estimated from a fit of curves of constant velocities to the 

H I velocity field 



s . We assume / to be unity. T^pm is an unknown quantity and we adopt a value of 100 K, typical of cold clouds in our Galaxy. 
The Hi column densities are then 1.5(±0. 2) x 10^' cm and 6.0(+1.5) x 10^*^ cm for the two components respectively. 



Table 4. H i and OH absorption results 



Parameter 


Hi 


OH(1667) 


OH(1665) 


Tl 


0.073 ± 0.004 


0.023 ± 0.003 


0.014 ± 0.004 


I'l (km s"') 


4705 ± 5 


4721 ±6 


4710 ± 19 


cr,,, (kms"') 


44.2 ± 4.4 


44.2 


44.2 




0.032 ± 0.004 


0.011 ±0.003 


0.006 ± 0.004 


V2 (km s"') 


4579 ± 10 


4585 ± 12 


4601 ± 45 


cr,,, (km s"') 


43.4 ±9.1 


43.4 


43.4 



4.4. OH absorption 

The OH spectrum shown in Fig. ^ was extracted at the radio position of Mrk 1 from the image cube made using the VLA data 
with a resolution of 5'.'43 x 4'.'93 x 27 km s The cube has an rms of 0.5 mJy beam 
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Fig. 7. The GMRT spectrum showing H i absorption in Mrk 1 . The dotted curve is the Gaussian fit to the absorption spectrum. 
The fitted parameters are given in Table The vertical lines along the velocity axis mark the positions of several velocity systems 
as indicated. 
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Fig. 8. The VLA spectrum showing OH absorption from Mrk 1. The velocity axis corresponds to the 1667 MHz line. The dotted 
curve shows the model spectrum of the 1665 and 1667 MHz OH lines. The fitted parameters are given in Table 0. In this velocity 
definition, the 1665 MHz line will appear at +360 km s ' from the 1667 MHz line. Two vertical lines in the velocity axis mark 
the expected positions of 1665 MHz absorption corresponding to the 1667 MHz absorption detected near 4721 and 4585 km s 



The velocity axis of Fig. || corresponds to the 1667 MHz OH line. In this velocity system, the 1665 MHz line will appear 
at +360 km s ' from the 1667 MHz line. Since the spectrum of Fig. ^ does not have enough baseline for the 1667 MHz line 
and not enough signal to noise ratio for both the 1665 MHz line to get reliable estimates for velocity dispersions of individual 
components, only peak optical depths and center velocities were fitted while the velocity dispersions were fixed at those values 
found in fitting the Hi absorption profile. This is a reasonably good assumption since both the 1665 and 1667 MHz profiles 
are qualitatively similar to the H i absorption profile. This procedure gave a reasonably good, though not unique, fit to the OH 
spectrum. The fitted parameters are given in Table 

The peak optical depths of the two components of the 1667 MHz absorption are 0.023+0.003 and 0.01 1 ±0.003 respectively. 
The column density of OH is estimated using the relation A^OH = 2.35 x 10'^ x {T^Jf) / T1667 dv cm where Tex is excitation 
temperature which is assumed to be 10 K and J Tig^v dv is the velocity integrated optical depth of the 1667 MHz line in units of 
km s"'. The OH column densities are estimated to be 6.0(+1.0) x 10'^ cm"^ and 2.9(±0.9) x lO'*" cm"^. The peak optical depth 
ratio ti667/ti665 of the stronger OH component is 1 .6 ± 0.5, indicating that this component is excited under LTE conditions - the 
ratio is predicted to be in between 1.0 and 1.8 for LTE excitations. This ratio for the weaker component, viz., 1.8 + 1.3 indicates 
that this is also most likely thermally excited. 
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Table 5. Interaction properties of the Mrk 1-NGC 451 system 



Parameter 


Value 


Projected velocity difference (km s"') 


117 


Projected separation (kpc) 


38 


Total dynamical mass (Mg) 


~ 10" 


Tidal radius (Mrk 1) (kpc) 


24 


Tidal radius (NGC 451) (kpc) 


29 


Impact parameter 


-0.1 


Dynamical friction time (Gyr) 


-0.2 


Orbital time (Gyr) 


-2.0 



5. Discussion 

5.1. Interaction of Mrk 1 with NGC 451 



The disturbed Hi morphology of Mrk 1 (Fig |]) indicates a gravitational interacti on possibly with the neare st companion 



NGC 45 1 . We explore this possibility using the two body interaction described in [Binney & Tremaine (1987)[ The dynami- 
cal masses of galaxies are estimated using rotation curves. Since the H i morphology of Mrk 1 is disturbed, it was not possible to 
obtain a reliable H i rotation curve. We used the Ha rotation curve of Mrk 1 (Keel 1996). The dynamical mass of NGC 45 1 was 
estimated using the H i rotation curve. The dynamical masses of Mrk 1 and NGC 451 are 3.4 x 10'" Mq and 4.5 x 10'" Mq 
respectively. The interaction parameters are listed in Table |[ 

The projected velocity difference between Mrk 1 and NGC 45 1 of ~ 117 km s indicates a minimum dynamical mass of this 
pair to be ~ 10" M©. This value of dynamical mass is in close agreement with the dynamical masses of Mrk 1 and NGC 451, 
indicating that Mrk 1 and NGC 451 are most likely in a bound system. Tidal radii (cf. eq. 7-84, Binney & Tremaine 1987) for 
the given masses of Mrk 1 and NGC 45 1 indicate that the outer regions of the H i disk of Mrk 1 can be perturbed easily. The 
dynamical friction time scale (cf. eq. 7-26 and 7-13b, Binney & Tremaine 1987) of - 0.2 Gyr for this system is much smaller 
than their orbital time scale of ~ 2 Gyr. This implies that the interaction is bound and will lead to a merger within a small fraction 
of the orbital time period of the two galaxies. 



5.2. Comparison of H i and OH velocities with other velocity systems 



Comparison with optical line velocities- From the comparison of the Hi and OH absorption velocities of Mrk 1 with the optical 
line velocities listed in Table |l], it appears that the higher velocity components of the H i and OH absorption are consistent with 
the [Oii]/l3727 line velocity, but are blueshifted by ~ 100 km s"' from the [Oiii]/15007, [Oi]/I6300 and the systemic velocity. 
We explain this discrepancy in terms of co-existence of photo-ionized and shock ionized gas in active galaxies. The [Om]/15007 
line is primarily due to excitation from a hard continuum, and therefore, should be arising close to the nucleus. The [Oii]/l3727 
line intensity is enhanced in shock ionized regions (Dopita & Sutherland 1995). Most often, optical line profiles are asymmetric 
and only peak line velocities are quoted without fitting a line profile. Mrk 1 is known to be such a case (Bergeron & Durret 1987, 
Dickinson et al. 1976). Such an analysis of the optical spectrum may bias the line velocities of different species toward different 
regions, e.g., the peak of the [On] line may indicate a region which is shock ionized while the [Oiii] line velocity may indicate 
gas which is close to the nucleus. We speculate that the higher velocity Hi and OH absorption component in Mrk 1 arises in a 
region which is pushed forward due to shocks, thereby giving an apparent blueshift to H i, OH, and [O ii] lines. The fact that the 
[O i] line velocity is close to the [O iii] line velocity, and hence associated with photo-ionized regions, is not surprising since the 
[Oi] line intensity is suppressed in the shock excited regions (Dopita & Sutherland 1995). 

Comparison with H2O megamaser and CO emission- The water megamasers are seen from Mrk 1 at a velocity of 4868 km s~' 
(Braatz et al. 1994). Since these masers are redshifted by ~ 90 km s"' from the systemic velocity of the galaxy, they are most 
likely the high velocity "satellite" features commonly seen in water megamaser galaxies (Braatz etal. 1997) and thought to 
originate in the accretion disks near the nuclei (Neufeldet al. 1994i. Since the observed Hi and OH absorption velocities in 
Mrk 1 are blueshifted from both the water megamaser velocity and from the systemic velocity of the galaxy, the absorption in 
the present case is most likely not related either to the g as in the accretion disk or to the torus close to the nucleus. The mean 
velocity of the CO emission from Mrk 1 is 4850 km s ' ( Vila-Vilaro et al. 1998 ), which is redshifted by ~ 150 km s ' from Hi 
and OH absorption velocities, implying that the gas traced via H i and OH absorption in Mrk 1 is also not related to the molecular 
gas traced by CO emission. 
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Table 6. Properties of the absorbing gas 







ly [\n If v^uiii J 






ft n(+ 1 'S'l X 1 
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7 sc+1 n't V 




9 1 c+o (S'l V 1 n^' 




O.Ot^XJ.^^ A 1\J 


^Hl/T^OH, 1667 


7 9 -I- n A 


Thi/Tqh, 1667 K^) 




' UH, loo// ' UH, loo3 


1.6 ± 0.5 


"^OH, 1667 /TqH, 1665 (2) 


1.8 ± 1.3 


Mo,al (H2)'t (cm-2) 


~ 1023 


* total V.*-'^ ^// ^ * total ^ ^/ 


~ 4.3 X 10"'' 


Av (mag) 


9-90 


A'h cm"^ 


~ 1.1 X 10" 


^IkcV 
photoelectric 


-30 


t : Assuming rspin = 100 K; X 
Assuming OH/H2 = 10-'' 


Assuming ri;x=10 K; tt : 



5.3. Kinematics and composition of tlie absorbing gas 

The general properties of the gas seen in absorption are summarized in Table |[ The total column density of OH is comparable 
to that observed in other active galaxies (e.g., Schmelz et al. 1986, Baan et al. 1985, 1992). Both components of the OH tran- 
sitions appear to be thermally excited as their optical depth ratios (t"i667/ti665) are between 1.0 and 1.8; the values predicted 
for excitations in LTE conditions. The optical depth ratios thi/tqh for both the components of the absorbing gas are ~ 3. This 
ratio has been found to be varying from as low as 5 to more than 400 in Galactic clouds (Dickey et al. 1981). The smaller values 
correspond to the molecular clouds while larger values correspond to the diffuse clouds. It is therefore suggested that the H i and 
OH absorption, in the present case, are associated with dense molecular clouds. 

The observed velocity dispersion (cr) of the 1667 MHz OH absorption, viz., 44 km s ' is higher than the typical velocity 
dispersions (cr = 3-7 km s ') in giant molecular clouds (GMCs) of the Galactic disk. However, several high dispersion clouds 
(cr ~ 40 km s ') have been detected in 18 cm OH absorption within a kpc of the Galactic centre (Boyce & Cohen 1994). The 
simplest explanation for such a high velocity dispersion could be a chance alignment of several normal GMCs along the line 
of sight, but Kumar & Riffert (1997)| have shown that the probability of such alignments is small. Alternatively, if the velocity 
dispersion is due to a single gravitationally bound system in virial equilibrium, the mass of such an object (assuming a size of 
50 pc) could be ~ 10^ Mq. Cloud-cloud collisions ( |Klein et al. 1994b| ) and interaction of shock with ISM ( |Klein et al. 1994a| ) 
are also known to enhance the internal velocity dispersions of molecular clouds. 

The OH column density is known to correlate with the visual extinction. Ay, of molecular clouds in our Galaxy 
( [Magnani et al. 1988| ). Magnani et al. (1988) found that N(OH)/Av ratios are in the range of 10'"* - 10'^ cm"^ mag"'. For the OH 
column density toward Mrk 1, these ratios indicate 9 < Ay < 90 toward the line of sight of Mrk 1. In comparison, Veilleux et al. 
(1997), based on some infrared measurements, obtained a lower limit on Ay to be 26 consistent with above predictions. 

Using OH/H2 = 10"^ (Liszt & Lucas 1996), the implied column density of H2 is ~ 10^^ cm"^. Using values of the photo- 
electric absorption cross sections from fVIorrison & McCammon (1983) for a gas having the solar abundance, a total hydrogen 
column density of ~ lO^^* cm"^ indicates that the optical depth for X-ray absorption at 1 keV will be ~ 30. Such a high value of 
the optical depth will absorb almost all soft X-radiation from the nucleus of Mrk 1. Consistent with this prediction, Mrk 1 has 

10^1 ■ 



not been detected as a X-ray source down to a sensitivity of - 



ergs (Fabbiano et al. 1992). 



6. Conclusions 

We have presented the observations of the Seyfert 2 galaxy Mrk 1 in the Hi 21 cm line using the GMRT and in the OH 18 cm 
line using the VLA. Unlike the optical morphology, the H i emission morphology of Mrk 1 indicates that this galaxy is disturbed 
which we interpret as due to tidal interactions with the nearby companion NGC 45 1 . We also showed based on the dynamical 
study of Mrk 1 - NGC 45 1 system that the interaction is bound leading to a merger within a small fraction of their orbital time 
period. This is consistent with the hypothesis that the bound interactions should be more efficient in triggering nuclear activities 
than unbound interactions. The H i and OH absorption detected toward the nucleus of Mrk 1 indicates an outflow of both atomic 
and molecular gas. The column densities of the detected H i and OH absorption indicate that the line of sight toward the nucleus 
of Mrk 1 is rich in both atomic and molecular gas. The gas detected in absorption is kinematically different than that traced via 
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CO emission and water megamaser emission from Mrk 1. We found evidences that shocks (presumably due to nuclear activities) 
can affect the kinematics of gas near the nucleus. The H i and OH absorption being blueshifted from the systemic velocity and the 
[Om]/I5007 velocity while consistent with the [O n]A3727 velocity is understood in terms of the shock ionization of gas (which 
predicts enhancement of the [O ii] line intensity) and an outflow of ISM in front of the shock. Based on the optical depth ratios 
and the line widths of the H i and OH absorption, we speculate that the absorption is arising in turbulent molecular clouds of 
similar types as those found near the Galactic centre. These observations also imply that the Une of sight toward the nucleus of 
Mrk 1 is heavily obscured. 
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